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This  report  looks  at  the  interaction  of  radiated  electromagnetic  fields  with 
earth  ground  in  military  or  law-enforcement  applications  of  high-power 
microwave  (HPM)  systems.  For  such  systems  to  be  effective,  the  microwave 
power  density  on  target  must  be  maximized.  The  destructive  and  construc¬ 
tive  scattering  of  the  fields  as  they  propagate  to  the  target  will  determine 
the  power  density  at  the  target  for  a  given  source.  The  question  of  field 
polarization  arises  in  designing  an  antenna  for  an  HPM  system.  Should 
the  transmitting  antenna  produce  vertically,  horizontally,  or  circularly  po¬ 
larized  fields?  Which  polarization  maximizes  the  power  density  on  target? 
This  report  provides  a  partial  answer  to  these  questions.  The  problems  of 
calculating  the  reflection  of  uniform  plane  wave  fields  from  a  homogeneous 
boundary  and  calculating  the  fields  from  a  finite  source  local  to  a  perfectly 
conducting  boundary  are  relatively  straightforward.  However,  when  the 
source  is  local  to  a  general  homogeneous  plane  boundary,  the  solution  can¬ 
not  be  expressed  in  closed  form.  An  approximation  usually  of  the  form 
of  an  asymptotic  expansion  results.  Calculations  of  the  fields  are  provided 
for  various  source  and  target  locations  for  the  frequencies  of  interest.  The 
conclusion  is  drawn  that  the  resultant  vertical  field  from  an  appropriately 
oriented  source  antenna  located  near  and  above  the  ground  can  be  signif¬ 
icantly  larger  than  a  horizontally  polarized  field  radiated  from  the  same 
location  at  a  1.3  GHz  frequency  at  observer  locations  near  and  above  the 
ground. 
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1.  Introduction 


Effective  military  or  law-enforcement  applications  of  high-power  mi¬ 
crowave  (HPM)  systems  in  which  the  HPM  system  and  the  target  system 
are  on  or  near  the  ground  or  water  require  that  the  microwave  power  den¬ 
sity  on  target  be  maximized.  The  power  density  at  the  target  for  a  given 
source  will  depend  on  the  destructive  and  constructive  scattering  of  the 
fields  as  they  propagate  to  the  target.  Antenna  design  for  an  HPM  sys¬ 
tem  includes  addressing  the  following  questions  about  field  polarization: 
Should  the  fields  the  transmitting  antenna  produces  be  vertically,  hori¬ 
zontally,  or  circularly  polarized?  Which  polarization  maximizes  the  power 
density  on  target?  (The  question  of  which  polarization  best  couples  to  the 
target  is  beyond  the  scope  of  this  report.)  While  this  report  does  not  com¬ 
pletely  answer  these  questions,  it  addresses  the  interaction  of  the  radiated 
electromagnetic  fields  with  earth  ground.  It  is  assumed  that  the  transmit¬ 
ting  antenna  and  the  target  (or  receiver)  are  located  above,  but  near  the 
surface  of  a  flat  idealized  earth  (constant  permittivity,  e,  and  conductivity, 
a)  ground.  First  an  ideal  vertical  dipole  (oriented  along  the  2-axis  perpen¬ 
dicular  to  the  ground  plane)  is  addressed.  The  horizontal  dipole  (parallel 
to  the  ground  plane)  follows. 
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2.  Problem  Formulation 


The  problems  of  calculating  the  reflection  of  uniform  plane  wave  fields 
from  a  homogeneous  boundary  and  calculating  the  fields  from  a  finite 
source  local  to  a  perfectly  conducting  boundary  are  relatively  straightfor¬ 
ward.  However,  when  the  source  is  local  to  a  general  homogeneous  plane 
boundary,  it  is  found  that  the  solution  cannot  be  expressed  in  closed  form. 
An  approximation  usually  of  the  form  of  an  asymptotic  expansion  results. 
The  problem  of  an  ideal  dipole  over  a  homogeneous  half-space  has  been  the 
topic  of  a  number  of  studies  starting  at  the  turn  of  the  last  century  with  the 
solution  provided  by  Sommerfeld  (1949)  and  leading  to  the  more  contem¬ 
porary  work  of  Banos  (1966)  and  King  et  al  (1994, 1992).  References  such 
as  Maclean  and  Wu  (1993)  address  in  detail  the  many  approaches  to  solv¬ 
ing  the  problem.  Considerable  controversy  has  surrounded  these  studies. 
We  will  not  attempt  to  derive  the  solution  or  otherwise  discuss  the  solu¬ 
tion  of  the  problem  in  this  report.  We  will  rely  on  the  work  of  King  et  al 
for  a  complete  and  concise  formulation  of  the  problem.  Figure  1  depicts  the 
problem  geometry.  The  King  expressions  have  been  encoded  and  solved  in 
MATLAB®  (1984-1999).  Comparisons  of  the  field  structures  are  provided 
for  various  source  and  target  locations  for  frequencies  of  interest. 

The  expressions  that  follow  are  constrained  by  the  magnitude  of  the  wave 
numbers  ( k  =  u^/xe)1/2,  u;  =  2i rx  frequency)  in  each  region: 

I&1.I  >  3  I&2I  >  (1) 

where,  after  King,  the  subscript  2  denotes  the  upper  half-space  (air)  and  the 
subscript  1  denotes  the  lower  half-space  (earth).  Also,  ft  is  the  permeability 
in  henries  per  meter,  e  is  the  permittivity  in  farads  per  meter,  and  a  is  the 
conductivity  in  siemens  per  meter.  The  subscripts  0  and  r  represent  free 
space  and  relative  to  free  space,  respectively. 
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2.1  Vertical  Dipole  Over  Earth 

The  electromagnetic  fields  from  a  dipole  with  dipole  moment  It  oriented 
perpendicular  to  and  at  a  height  z  —  d  above  the  ground  plane  have  three 
components  in  cylindrical  coordinates.  The  magnetic  field  is  symmetric 
about  the  2-axis,  perpendicular  to  the  direction  of  propagation.  For  that 
reason  the  fields  will  be  termed  transverse  magnetic  (TM)  fields.  (We  will 
find  that  the  horizontal  dipole  has  TM  and  transverse  electric  (TE)  com¬ 
ponents.)  The  formulation  as  provided  in  King  et  al  (1994)  is,  referring  to 
figure  1, 


H(j){p,  z)  — 


Ep(p,  z) 


2*  (ifc)  *«-"■*■(*>) 

'  eifc2ri  ( _£_\  ( z=d\  (ikz  _  3 _ 3i\ 

2  \rij  \  n  J  \  n  rf  k2r\ ) 

Ujpolt  ,eik2r2  (  n  \  ( z+d\  (ikz  3  3 i  \ 

\  +—  [r2)  {—)  ^  -  to* ) 

r ( _£_\  (ik 2 _ ±-\  _  ( (  2L  ^ 
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~k2X I /  V  r2  /  V  r2  rf  k2r%  J 

^fe-iPF(P) 


(2) 

(3) 

(4) 


where 

n  =  [p2  +  {z-  d)2] 2  , 

—  [p2  +  (^  +  d)2  5  ,  (5) 

p  =  fc2r2  fk2r2  +  k  1  (z  +  d)  \  2  ^ 

V  &2/>  / 

POO  pit  1 

F(P)=  -—-rdt=-(l  +  i)-C2(P)-iS2(P),  (7) 

dp  (27rt)  2  2 

and  C2(P),S2(P)  are  the  Fresnel  integrals  as  defined  in  Abramowitz  and 
Stegun  (1970).  Since 


we  have 


C2{P)  +  iS2{P)=l-{l+i)~ 


where  w(z),  called  the  plasma  dispersion  function  with  complex  argument, 
is  a  form  of  the  error  function  defined  in  Abramowitz  and  Stegun  (1970).  A 
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convenient  MATLAB®  m-file  is  available  for  solving  the  function  with  com¬ 
plex  argument  (Chase)  and  is  provided  in  the  appendix.  Note  that  \Ti  can 
be  written  as  and  that  F(P)  always  appears  with  e~iP ,  thus  cancel¬ 
ing  the  exponential  term  in  equation  (9).  This  is  reflected  in  the  MATLAB® 
m-files  that  calculate  the  fields  (provided  in  the  appendix). 

The  equations  have  been  written  such  that  the  direct  and  reflected  compo¬ 
nents  appear  first.  The  last  term  is  referred  to  as  the  surface  or  lateral  wave. 
Often  it  is  called  the  Norton  surface  wave  from  the  engineering  models  he 
developed  in  the  mid-1950s.  The  equations  reduce  to  the  fields  above  a  per¬ 
fectly  conducting  ground  when  k\  —  •  oo;  the  surface  or  lateral  wave  term 
then  goes  to  zero. 

2.2  Horizontal  Dipole  Over  Earth 

The  electromagnetic  fields  produced  by  an  ideal  dipole  oriented  at  a  height 
z  —  d  above  and  parallel  to  the  ground  plane  consist  in  general  of  both  TM 
and  TE  components.  The  formulation  as  provided  in  King  et  al  (1992)  for 
the  TM  wave  components  in  cylindrical  coordinates,  referring  to  figure  1 
and  the  above  definitions,  is 


The  TM  fields  are  zero  broadside  to  the  dipole  orientation,  <f>  =  f .  The  TE 
components  are 


4 


Hp{p,  <f>,  z) 


Hz(p,  4>,  z) 


E<p(p ,  4>,  z) 


These  are  the  predominant  fields  broadside  to  the  dipole  orientation.  Here 
Etj,  is  often  termed  the  horizontal  electric  field.  Again,  the  equations  reduce 
to  the  fields  above  a  perfectly  conducting  ground  when  k\  — >  oo;  the  surface 
or  lateral  wave  term  goes  to  zero. 


2.3  Comparisons 


The  peak  power  radiated  by  a  unit  dipole  is  given  in  Collin  and  Zucker 
(1969): 


p  _  fcjCo 
127 r  ’ 


(16) 


where  Co  is  the  free-space  impedance  (y/^)  •  The  fields  and  power  density 
comparisons  that  follow  are  normalized  to  one  watt  radiated  peak  power. 
To  obtain  the  unnormalized  quantities,  simply  multiply  the  power  results 
by  P  and  the  field  results  by  Pz.  The  outward  component  of  the  complex 
Poynting  vector  over  a  closed  surface  is 


V2  —  *  —  ‘  —  ^complex  )  (17) 

S 


where  the  negative  sign  indicates  power  flow  away  from  the  surface.  Our 
interest  is  in  the  real  part  of  the  power  density  at  the  observer  (or  target) 
location.  In  our  cylindrical  coordinate  system  for  the  TM  components,  this 
becomes 


1/2  Re  (E  x  H*)  =  1/2  Re  (p0EzH^  *  +  zqE^  *)  ,  (18) 
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where  p0  and  z0  are  the  unit  vectors  in  cylindrical  coordinates.  Near  the 
ground,  the  power  flow  is  predominantly  radial  with  a  small  2  component. 
And,  for  the  TE  components,  we  have 

V2  ^  (E  x  H*)  =  1/2  Re  (p()E4>Hz  *  -  z0E4>Hp  *)  .  (19) 

For  the  vertical  dipole,  the  power  density  at  the  observer  (on  target)  will  be 
Pv  =  1/2  (Re  (E$HZ  *)  +  Re  (E^H,  *))  .  (20) 

The  horizontal  dipole  in  general  will  produce  a  target  power  density  of 
Ph  =  V2  ({Re  (E^H,  *)  -  Re  (ESH+  *)}  +  {Re  (. E pH^  *)  -  Re  (E+Hp  *)})  .  (21) 

For  broadside  calculations  this  becomes 

Ph  =  V2  (Re  (E<t>Hz  *)  +  Re  (E^HP  *))  .  (22) 

The  calculations  that  follow  are  limited  to  a  frequency  of  1.3  GHz  (wave¬ 
length,  Ao,  is  0.23  m)  and  dipole  heights  of  1  to  3  m.  Observer  (target) 
heights  range  from  0  to  5  m.  The  frequency  of  1.3  GHz  is  chosen,  since  most 
of  the  HPM  source  and  antenna  design  work  at  ARL  is  centered  around  that 
frequency.  The  height  ranges  are  chosen  to  be  consistent  with  ground  vehi¬ 
cle  source  and  target  applications.  Five  classes  of  ground  parameters  that 
are  representative  of  distinctly  different  terrain  will  be  addressed.  These 
five  classes  are  those  discussed  in  King  et  al  (1994)  and  are  given  in  table  1. 


2.3.1  Comparison  of  Field  Components  Near  the  Earth 


Table  1.  Earth 
parameters. 

When  thinking  about  the  interaction  of  electromagnetic  waves  with  the 
earth's  surface,  we  often  view  the  earth  as  a  perfect  conductor.  The  hori¬ 
zontally  polarized  field  is  reflected  with  the  opposite  sign  of  the  incident 
field,  leading  to  a  difference  or  destructively  interfered-with  field.  The  ver¬ 
tically  polarized  field  is  reflected  with  the  same  sign  as  the  incident  field, 
leading  to  a  sum  or  constructively  interfered-with  field.  Figure  2  compares 
the  resultant  primary  field  components  produced  from  a  horizontal  and  a 
vertical  ideal  dipole  located  at  the  same  point  in  space  over  a  perfectly  con¬ 
ducting  ground  plane.  From  these  calculations,  one  might  conclude  that 

Case  cr*  S/m  er* 

1  Sea  water  4  80 

2  Wet  earth  A  12 

3  Dry  earth  .04  8 

4  Lake  water  .004  80 

5  Dry  sand  .000  2 

*The  variable  name  representing  a 
in  the  MATLAB  m-files  is  SIGMA 
and  for  er  it  is  EPSREL.  These  vari¬ 

able  names  appear  on  many  of  the 
figures. 
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Figure  2.  Comparison  of 
main  electric  field 
components  from  a 
vertical  (Ez)  and  from  a 
horizontal  ( E $) 
broadside  ideal  dipole 
over  a  perfectly 
conducting  ground.  Ez 
is  a  constructively 
interfered-with  field, 
while  E<p  is  a 
destructively 
interfered-with  field. 


Normalized  field  magnitude  (V/m) 


an  antenna  that  radiated  fields  polarized  such  that  the  electric  field  was 
essentially  perpendicular  to  the  ground  plane  would  produce  the  largest 
fields  and,  consequently,  the  greater  power  densities  on  a  target  or  near 
the  ground.  This  is  generally  not  the  case  for  real  earth  grounds.  While 
the  perfectly  conducting  ground  plane  model  may  provide  some  insight 
into  the  interaction  with  horizontally  polarized  fields,  it  clearly  is  an  inad¬ 
equate  model  for  vertical  polarization.  The  complex  reflection  coefficients 
for  a  plane  wave  polarized  with  the  electric  field  in  the  plane  of  incidence 
(vertical  polarization)  and  with  a  plane  wave  polarized  with  the  electric 
field  perpendicular  to  the  plane  of  incidence  (horizontal  polarization)  are 
given  by  the  following  Fresnel  expressions.  For  vertical  polarization,  it  is 
(Collin,  1985) 


Pcomplex  — 


{£r  li^o) 

1  sin  —  y 

|  —  cos2  ip 

|  sin  ip  +  y 

1  —  COS2  tp 

(23) 


and  for  horizontal  polarization,  it  is 


Pcomplex  — 


sin  'ip  —  ^ 

—  COS2  Ip 

sin  ip  +  ^ 

—  cos2  ip 

(24) 


where  ip  is  the  angle  of  incidence  as  measured  between  the  ray  path  and 
the  surface  of  the  earth  (see  fig.  1).  These  reflection  coefficients  for  the  five 
cases  of  table  1  are  presented  in  figure  3.  One  can  clearly  see  that  for  shal¬ 
low  grazing  angles,  a  Brewster  angle  effect  exists  (often  called  a  pseudo- 
Brewster  angle).  The  result  is  destructive  interference  of  the  vertically  po¬ 
larized  field  for  incident  angles  that  are  less  than  this  pseudo-Brewster 
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Figure  3.  Reflection 
coefficients  for 
vertically  and 
horizontally  polarized 
plane  waves  of  1.3-GHz 
frequency  incident  on  a 
flat  ground  for  five 
earth-parameter  cases 
listed  in  table  1.  Note 
that  case  5  provides  a 
true  Brewster  angle 
(since  conductivity  is 
zero)  and  case  4  is  very 
near  to  producing  a  true 
Brewster  angle  at 
1.3-GHz  frequency: 

(a)  vertical  polarization 
and  (b)  horizontal 
polarization. 


Incidence  angle  (°) 


angle.  This  would  lead  one  to  conclude  that  horizontally  and  vertically  po- 
larized  field  levels  above  a  ground  would  be  comparable  at  shallow  angles 
of  incidence.  The  Fresnel  equations  (22)  and  (23)  do  not  take  into  account 
the  surface  or  lateral  wave  produced  by  finite  sources  above  the  ground. 
They  are  also  not  applicable  to  near-field  problems.  We  can  compare  the 
reflection  coefficients  as  determined  from  the  Fresnel  equations  with  the 
results  of  our  complete  solution.  To  do  this  with  the  complete  equations 
(eqs  (2)  to  (4)  for  the  vertical  dipole  and  eqs  (10)  to  (15)  for  the  horizontal 
dipole),  we  fix  the  radial  (or  range)  distance  to  1000  m.  The  dipole  height  is 
varied  to  provide  a  range  of  incident  angles.  The  observer  will  be  fixed  at  a 
1-m  height.  The  incident  and  "reflected"  fields  are 

F'vtotai  =  +  E%  ,  (25) 

for  the  vertical  dipole  and 

•Ehtotal  =  \/-^J  +  +  El  ,  (26) 

for  the  horizontal  dipole.  For  broadside  calculations,  this  simply  becomes 
E<l>-  The  terms  involving  n  in  equations  (2)  to  (4)  and  (10)  to  (15)  repre¬ 
sent  the  incident  field  and  the  remaining  terms  represent  the  reflected  and 
surface  wave  terms.  The  reflection  coefficient  will  be  calculated  as  the  to¬ 
tal  incident  electric  field  divided  into  the  total  reflected  and  bound  (sur¬ 
face  wave)  fields.  Figures  4  and  5  provide  the  comparison  of  the  magni¬ 
tude  of  the  reflection  coefficient  for  the  five  cases  of  table  1.  For  the  most 
part,  the  Fresnel  expressions  are  a  good  approximation  for  the  calculation 
of  horizontally  polarized  field  values  (fig.  5)  above  realistic  earth  ground  at 
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Figure  4.  Fresnel 
reflection  coefficient  for 
a  vertically  polarized 
plane  wave  field 
compared  to  "reflection 
coefficient"  as 
calculated  by  complete 
formulation  for  a 
vertical  dipole  over 
homogeneous  earth  as 
observed  at  1-m  height 
1000  m  down  range  at  a 
frequency  of  1.3  GHz: 
(a)  sea  water,  (b)  wet 
earth,  (c)  dry  earth, 

(d)  lake  water,  and 

(e)  dry  sand. 
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1.3  GHz.  The  horizontally  polarized  fields  diverge  from  the  Fresnel  expres¬ 
sions  for  low  conductivities  and  relative  dielectric  constants.  The  resulting 
total  fields  in  such  cases  predicted  by  the  Fresnel  expressions  will  be  larger 
than  the  complete  solution  prediction.  The  results  for  vertical  polarization 
(fig.  4)  show  significant  divergence  for  shallow  incident  angles.  In  this  case 
the  resulting  total  fields  predicted  by  the  complete  solution  can  be  signifi¬ 
cantly  higher  than  that  predicted  from  employing  the  Fresnel  expressions. 
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Figure  5.  Fresnel 
reflection  coefficient  for 
a  horizontally  polarized 
plane  wave  field 
compared  to  "reflection 
coefficient"  as 
calculated  by  complete 
formulation  for  a 
horizontal  dipole  over 
homogeneous  earth  as 
observed  at  1-m  height, 
broadside  1000  m  down 
range  at  a  frequency  of 
1.3  GHz:  (a)  sea  water, 
(b)  wet  earth,  (c)  dry 
earth,  (d)  lake  water, 
and  (e)  dry  sand. 
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Figures  6  and  7  are  plots  of  the  fields  produced  by  the  two-dipole  orien¬ 
tations;  these  plots  compare  the  predominant  field  components  for  each 
dipole  orientation  at  1000  and  100  m,  respectively  The  plots  for  the  100-m 
case  (fig.  7)  show  the  beginning  of  the  lobe  effect  produced  by  the  phase 
difference  between  the  incident  and  reflected  wave.  This  is  shown  more 
dramatically  in  figure  8. 


Figure  6.  Comparison  of 
principal  fields  from  an 
ideal  dipole  oriented 
perpendicular  and 
horizontal  to  a 
homogeneous  flat  earth. 
In  each  case,  dipole  is 
placed  2  m  above 
ground  plane  and 
observer  or  target  is 
1000  m  down  range: 

(a)  sea  water,  (b)  wet 
earth,  (c)  dry  earth, 

(d)  lake  water,  and 

(e)  dry  sand. 
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Figure  7.  Comparison  of  (a) 
principal  fields  from  an 
ideal  dipole  oriented 
perpendicular  and  1 

horizontal  to  a  S 

0) 

homogeneous  flat  earth,  t 

In  each  case,  dipole  is  j 

placed  2  m  above  ° 

ground  plane  and 
observer  or  target  is 
100  m  down  range: 

(a)  sea  water,  (b)  wet 
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Figure  8.  Effect  of 
ground  reflection  on 
primary  field 
components  near 
ground  for  typical  earth 
parameters.  Phase 
difference  between 
incident  and  reflected 
waves  results  in 
development  of  field 
lobes  that  are  more 
pronounced  as 
radiating  antenna  is 
approached:  (a)  100  m 
down  range  and  (b) 
1000  m  down  range. 
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3.  Conclusion 


A  suite  of  MATLAB0  m-files  have  been  developed  to  calculate  the  electro¬ 
magnetic  fields  produced  by  a  vertical  and  a  horizontal  infinitesimal  unit 
dipole  over  a  homogeneous  flat  (ground)  plane.  Calculations  for  the  fields 
above  the  ground  plane  have  been  made  for  various  ground-plane  conduc¬ 
tivities  and  relative  dielectric  constants.  The  calculations  bound  the  practi¬ 
cal  range  of  parameters  representative  of  natural  earth  terrain. 

For  a  frequency  of  1.3  GHz,  where  the  dipole  and  the  observer  are  close  to 
the  ground  plane  (<3  m),  significant  difference  is  seen  in  the  magnitude  of 
the  fields  from  either  dipole  orientation.  The  power  density  on  target  will 
be  much  larger  for  vertical  dipole  orientation.  The  effects  of  rough  terrain, 
foliage,  or  scattering  from  manmade  or  natural  objects  in  the  path  from  the 
dipole  to  target  may  alter  this  conclusion.  How  these  other  scatterers  might 
affect  the  field  structure  at  a  target  at  1.3  GHz  is  not  known  at  this  time.  If 
the  effects  are  random  in  nature,  the  present  conclusion  is  most  likely  still 
valid.  From  the  simple,  smooth,  flat  ground  model,  then,  one  must  con¬ 
clude  that  vertical  polarization  (antenna  radiating  a  vertically  polarized 
field)  will  deliver  the  most  energy  to  the  target.  Unless  the  target's  pref¬ 
erence  for  field  orientation  for  maximum  pickup  is  known,  a  vertically  po¬ 
larized  antenna  may  in  fact  be  the  best  choice  for  a  ground  weapon  system. 
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Appendix.  MATLAB®  m-Files 


This  appendix  documents  the  MATLAB®  m-files  that  implemented  the  field 
equations  for  the  vertical  and  horizontal  dipole  and  array  cases.  Slight  mod¬ 
ification  may  be  required  to  obtain  all  the  calculations  presented. 


WERF  Function 

function  w  =  werf(z,N) 

X  WERF(Z,N)  Plasma  Dispersion  Function  with  complex  argument. 

X 

X  Computes  the  function  w(z)=:exp(-z''2)*erf c(-iz)  using  a  rational 
X  series  with  N. terms.  N  should  be  a  power  of  2  or  it  gets  SLOW. 

X  Default  value  of  N  is  64.  z  can  be  a  matrix  of  values. 

*/,  Taken  from  Siam  Journal  on  Numerical  Analysis,  Oct  1994,  V31  ,#5 . 
*/,  Modified  by  R.  Chase  to  work  for  all  z  (4/3/95). 

X 

X  N=32  gives  approx  14  place  accuracy,  N=64  is  better  and 
X  it  seems  as  fast. 

X 

X  See  wfn.m  --  Draws  graph  in  Abramowitz  &  Stegun, 

X  Handbook  of  Mathematical  Functions,  p.  298 


if  nargin  ==  1,  N=64;  end 

X  Default  value  for  N 

M=2*R;  M2=2*M;  k=[-M+l:l:H-l] *; 

X  K2=no.  of  sampling  points 

L=sqrt(N/sqrt (2) ) ; 

X  Optimal  choice  of  L 

theta=k*pi/M;  t=L*tan(theta/2); 

X  Variables  thetaand  t 

f=exp(-t .  ~2) .  *(Ll+t .  ~2) ;  f  =  [0 ;f] ; 

X  function  to  be  transformed 

a=real(f f t (ff tshift (f ) ) ) /M2 ; 

X  Coefficients  of  transform 

a=flipud(a(2:N+l)); 

X  Reorder  coefficients 

nz=imag(z)  <=  0; 

X  Find  im(z)  <=0 

•z(nz)=conj(z(nz)); 

X  Use  conj  for  above 

Z=(L+i*z) . / (L-i*z) ;  p=polyval(a,Z) ; 

X  Polynomial  evaluation 

w=2*p . /(L-i*z)  /2+(l/sqrt (pi) ) . / (L-i*z) ; 

X  Evaluate  w(z) 

w(nz)=2*exp(-(conj(z(nz))/2))  -  conj(w(nz)); 

X  Handle  im(z)  <=  0 

Xif  all(imag(z)==0),  w=real(w);  end 

X  Rtn  real  if  real 
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Constants 


function  [w , cv, epso , uO , zo , kO , kl , kappap , lo] =cnstdg(f ) 

'/■  This  m-file  contains  the  basic  constants  to  be  used  by  various 
'/,  m-files  and  functions 

l 

'/,  f  is  input  in  GHz 

'/,  global  EPSREL  SIGMA 

w=2*pi*f*le9; 

cv=2.99792458e8; 

epso=(l/(43*:pi)  )*le7*(l/cv"2) ; 

epsl=EPSREL; 

u0=4*pi*le-7; 

zo=sqrt(uO/epso); 

kO=w*sqrt(epso*uO) ; 

epsr= (epsl-j *SIGMA . / (w*epso ) ) ; 

ki=kO*sqrt(epsr) ; 

kappap=-kO/sqrt(epsr+l); 

lo=cv/(f*le9); 


Vertical  Dipole 

function  [ez,er,hp]=dipg(f ,d,rho,z) 

l 

l  Electromagnetic  fields  from  a  vertical  z-directed  current  element 

7.  at  a  height,  d,  over  a  conducting  dielectric  plane  (ground)  calculated 

'/.  by  the  King/Sandler  model.  This  formulation  is  that  developed 

'/,  as  equations  6,  7,  and  8  of  "The  Electromagnetic  Field  of  a  Vertical  Electric 

'/,  Dipole  over  the  Earth  or  Sea",  IEEE  Trans,  on  Antennas  and  Propagation, 

'/.  March  1994,  Vol  42  No.  3,  page  383.  It  is  the  same  as  the  King/Owens/Wu 
'/,  formulation.  This  formulation  is  that  developed  as  equations  4.2.30-32  of 
'/,  "Lateral  Electromagnetic  Waves",  Springer-Verlag,  1992,  pp  293-297. 

I 

'/,  >  f-  frequency;  rho-  radial  distance  from  the  z-axis 

'/,  <  z-  height  of  observer  above  ground  plane  (input  as  an  array) 

'/■  >  ez,er,hp-  field  components  at  the  observer,  where  the  second  letter 

'/,  <  designates  the  component  —  z,  r  (rho),  p  (phi) 

l 

l  The  formulation  is  subject  to  Ikl(ground) |>3|k2(air) I 

l 

[w,c,eps0,u0,zo,k2,kl, kappap, lo]=cnstdg(f); 

l 


18 


'/,  King,  et.  al  assume  an  exp(-i*omega*time)dependency  thus  we  convert  kl 

7, 

kl=conj(kl); 

7,B2=(kl/k2)"2; 
k21=k2/kl ; 

7, 

7,  We  assume  a  unit  dipole 

7, 

il=l; 

rl=(rho>(z-d).‘2).\5; 

r2=(rho'2+(z*d)."2).".5; 

sd=rho./rl; 

sr=rho./r2; 

cd*(z-d)./rl; 

cr=(d*z)./r2; 

p=(r2*k2‘3/(2*kl“2))  ,*((k2*r2+kl*(z+d) ) . /(k2*rho) ) .  “2; 

7. 

7,  The  attenuation  function  less  the  exp(i*p)  term  -  related  to  the 
7,  error  function  is 

7. 

Ferf=((l+i)/4)*werf(sqrt(i*p)); 

7, 

7.  Develop  the  terms  for  the  field  expressions 

7, 

tl=w*u0*il/(2*pi*k2); 

tll=-il/(2*pi); 

t2=exp(i*k2*rl)/2; 

t3=(i*k2 .  /rl)-(l . /rl . • 2)-i . /(k2*rl .  ‘3) ; 
t31=(i*k2./rl)-(l./rl."2); 

t4=(cd.‘2).*((i*k2./rl)-(3./rl.‘2)-(3*i)./(k2*rl/3)); 
t41=((i*k2./rl)-(3./rl . “2)-(3*i) . /(k2*rl . ‘3)) ; 
t5=exp(i*k2*r2)/2; 

t6=(i*k2./r2)-(l./r2.‘2)-i./(k2*r2.‘3); 

t61=(i*k2./r2)-(l./r2.‘2); 

t7=(cr."2).*((i*k2./r2)-(3./r2."2)-(3*i)./(k2*r2."3)); 
t71=((i*k2./r2)-(3./r2.“2)-(3*i)./(k2*r2."3)); 
t8=(2*t5*(k2'3)/kl) . *((pi . / (k2*r2)) . * . 5) . *sr . *Ferf ; 
t81=((k2"3)/kl)*((pi .  /(k2*r2)) . " .  5) .  *Ferf ; 
t82=(2*t5.*(k2'3)/kl) . *((pi . /(k2*r2) ) . ‘ . 5) . *Ferf ; 

7, 

7.  Calculate  the  fields 
7,  ez=tl*(t2 .  *(t3-t4)+t5 .  *(t6-t7)-t8) ; 
er=-tl*(t2.*sd.*cd.*t41+t5.*sr.*cr.*t71-k21*2*t5.*(sr.*t61-t81)); 
hp=tll*(t2 . *sd . *t31+t5 . *sr . *t61-t82) ; 


Horizontal  Dipole 


function  [ez ,ep,er,hz,hp)hr]=dipgh(f ,d,phi ,rho,z) 

7, 

l  Electromagnetic  fields  from  a  horizontally  directed  (perpendicular  to  z) 

7.  current  element  located  in  the  phi=0  plane 

7,  at  a  height,  d,  over  a  conducting  dielectric  plane  (ground)  calculated 
l  by  the  King/Owens/Wu  formulation.  This  formulation  is  that  developed 
7,  as  equations  7.10.75,  7.10.80,  7.10.84,  7.10.92,  7.10.93  and  7.10.94  of 
7,  "Lateral  Electromagnetic  Waves",  Springer- Verlag,  1992,  pp  293-297. 

7, 

7,  >  f-  frequency(GHz);  phi-  angle  about  the  z-axis,  rho-  radial  distance  from 

7.  <  the  z-axis,  z-  height  of  observer  above  ground  plane  (input  as  an  array) 

7,  >  ez,ep,er,hz,hp,hr-  field  components  at  the  observer,  where  the  second 

7,  <  letter  designates  the  component  --  z,  r  (rho),  p  (phi) 

7, 

l 

l  The  formulation  is  subject  to  Ikl(ground) I >3|k2(air) | 

7, 

[w,c,eps0,u0;zo,k2,kl,kappap,lo]=cnstdg(f); 

7, 

7*  King,  et.  al  assume  an  exp(-i*omega*time)dependency  thus  we  convert  kl 

7, 

kl=conj (kl) ; 

7,N2=(kl/k2)A2; 

k21=k2/kl; 

k2p=k2/rho; 

l 

7,  We  assume  a  unit  dipole 

7, 

il=l; 

cp=cos(phi); 

sp=sin(phi); 

rl-(rhoA2+(z-d).A2).\5; 

r2=(rhoA2+(z+d).A2).\5; 

sd=rl/rho; 

sdo=l./sd; 

sr=r2/rho; 

sro=l./sr; 

cd=(z-d)./rl; 

cr=(d+z)./r2; 

p=(r2*k2A3/(2*kr2)).^((k2*r2+klt(z-d))./(k2+rho)).A2; 
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X 

7.  The  attenuation  function  less  the  exp(i*p)  term  -  related  to  the 

'/,  error  function  is 

'/,  Ferf=((ln)/4)*werf(sqrt(i*p)); 

'/, 

7,  Develop  the  terms  for  the  field  expressions 

X 

tl=»*u0*il*cp/(4*pi*k2); 
tll=-v*u0*il*sp/(4*pi*k2) ; 
tl2=il*sp/(4*pi); 
tl3=il*cp/(4*pi); 
t2=exp(i*k2*rl); 

t32=(i*k2 . /rl)-(l . /rl . ~2)-i . /(k2*rl . "3) ; 
t3=(2 . /rl . "2)+2*i ./(k2*rl . "3) ; 
t31=(i*k2./rl)-(l./rl."2); 

t4=(cd."2).*((i*k2./rl)-(3./rl.“2)-(3*i)./(k2*rl.'3)); 

t41=(cd).*((i*k2./rl)-(3./rl.‘2)-(3*i)./(k2*rl."3)); 

t5=exp(i*k2*r2); 

t62=(i*k2./r2)-(l./r2.‘2)-i./(k2*r2.‘3); 

t6=(2./r2.*2)*2n./(k2*r2.‘3); 

t61=(i*k2./r2)-(l./r2."2); 

t7=(cr."2).*((i*k2./r2)-(3./r2."2)-(3*i)./(k2*r2.“3)); 
t71=(cr).*((i*k2./r2)-(3./r2/2)-(3*i)./(k2*r2.-3)); 
t72=((l./r2."2M3*i)./(k2*r2."3)-3./((k2‘2)*r2.~4)); 
t73=(cr.A2).*((i*k2./r2)-(6./r2."2)-(15*i)./(k2*r2.“3)); 
t8=((k2"3)/kl)*C(pi • /(k2*r2)) .5) .*sr . *Ferf ; 
t81=((k2‘3)/kl)*((pi/(k2*r2)).".5).*Ferf; 
t82=(2*t5*(k2A3)/kl) . * ( (pi . /(k2*r2)) . * . 5) . *Ferf ; 
t83=((pi . / (k2*r2)) .  \  5) . *Ferf ; 

l 

l  Calculate  the  fields 

l 

.ez=tl*(-t2.*sdo.*t41+t5.*sro.*t71-2*k21*t5.*(sro.*t61-t81)); 
er=tl*(t2>(t3H4)-t5.*(t6H7-2*k21«6h2*(k2r2)*(t62-t8))); 
ep=tll*(t2.*t32-t5.n62-t5.*(-2*k21*cr.*t61+2*(k2r2)*(t6+t7). . . 
+2*i*k21"3*k2p*(sr. "2) .*t83)) ; 

hr=tl2*(t2.*cd.*t31-t5.*cr.*t61+2*k21*t5.*(t6+i*k21*k2p*(sr.A2).n83+t7)); 
hp=tl3*(t2.*cd.*t31-t5.*cr.*t61+2*k2i*t5.*(t62-k21*(k2‘2)*(sr).*t83)); 
hz=tl2*(t2 . *sdo . *t31-t5 . *sro . *t61*2*sro . *t5 . * (k21*t71 (k21“2)*(t72+t73)) ) ; 


Plot  m-File  for  Fields 

X 

X  This  m-file  plots  the  fields  over  a  conductive  flat  earth  produced  by  an  ideal 
/■  dipole  placed  a  distance  d  above  the  earth.  It  compares  the  results  from 
'/,  a  vertical  and  horizontal  dipole. 

X 

X 

X  Establish  the  problem  conditions 

x 

x 

'/.  EPSREL-  Relative  dielectric  constant;  SIGMA-  Earth  conductivity  (S/m) 

X 

EPSREL=80 ; SIGHA=4 ; 
global  EPSREL  SIGMA 
XEPSREL= 1 2 ; S IGMA= . 4 ; 

‘/.EPS  REL=8 ;  SI  GH  A= .  04 ; 

XEPSREL=80;SIGMA=.004; 

*/EPSREL=2 ;  SIGMA= .  000 ; 

X 

'/.  Location  of  dipole  (m) 

X 

d=2; 

y. 

X  Location  of  observer,  rho  (n);  phi  (radians);  z  (m)  — >  an  array 

y. 

rho=1000; 

phi=pi/2; 

z*.l*[l:l:30]; 

y. 

X  f-  frequency  in  GHz 

X 

f=1.3; 

X 

X  Field  normalization  factor  -  one  Watt  radiated 

y. 

[w,c,epsO,uO,zo,k2,kl,kappap,lo]=cnstdg(f); 

Fn=sqrt(12*pi/((k2"2)*zo)); 

X 

X  Horizontal  dipole  fields 
X 

[ez,ep,er,hz,hp,hr]=dipgh(f ,d,phi,rho,z); 

plot(Fn*abs(ep),z,’-b') 

hold 
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X 

X  Vertical  dipole  fields 

X 

[ezv , erv , hpv] =dipg(f , d , rho , z) ; 
plot(Fn*abs(ezv) ,z, 
ts=Fn*abs(erv(l)); 

title( ’Vertical  and  Horizontal  Ideal  Dipole  fields  over  Ground’) 
text(ts,z(19) , [’SIGMA  =  ’ ,num2str( SIGMA) , ’  S/m  EPSREL= 

’ ,num2str(EPSREL)] ) 

text(ts,z(18), [’Dipole  height  =  ’ ,num2str(d) , ’  meters;  Frequency  = 

’,num2str(f),’  GHz’]) 

text(ts,z(17), [’Range  =  ’,nuni2str(rho)}’  meters’]) 

ylabel( ’Observer  height  -  meters’) 

xlabel( ’Normalized  field  magnitude  -  V/m’) 

legend( ’Horizontal  E-field,  horizontal  dipole’, ’z-directed  E-field, 

vertical  dipole’) 

hold 

Fresnel  Reflection  Coefficients  and  Plots 

function  [rcomv , rcomh] =Fresnell (f ,psi ) 

X 

l  This  routine  calculates  the  Fresnel  reflection  coefficients  for 
l  vertical  and  horizontal  incident  fields 

X 

X  after  Collin,  "Antenna  and  Radiowave  Propagation",  page  345 

X 

X  f -frequency  in  GHz;  psi-  array  of  incident  angles 

X 

X  Code  returns  the  arrays  rcomv  and  rcomh,  the  vertical  and  horizontal 
X  complex  reflection  coefficients,  respectively. 

X 

global  EPSREL  SIGMA 
w=2*pi*f*le9; 
cv=2.99792458e8; 
epso=(l/(4*pi))*le7*(l/cv"2); 
epsr=(EPSREL-j*SIGMA . / (w*epso)) ; 

rcomv=(epsr*sin(psi)-sqrt(epsr-cos(psi)/2))./(epsr*sin(psi)+sqrt(epsr-cos(psi)/2)); 

rcomh=(sin(psi)-sqrt(epsr-cos(psi).~2))./(sin(psi)+sqrt(epsr-cos(psi)/2)); 


/•  This  m-file  plots  the  reflection  coefficient  for  vertical  and  horizontal 
l  fields  over  a  homogeneous  ground  plane 

l 

l 

I  EPSREL-  Relative  dielectric  constant;  SIGMA-  Earth  conductivity  (S/m) 

l 

EPSREL=80;SIGKA=4; 

f=1.3; 

global  EPSREL  SIGMA 

psi= [pi/1000 :pi/1000 :pi/2] ; 

trcomv,rcomh]=Fresnell(f ,psi) ; 

EPSREL=12;SIGMA=.4; 

[rcomvl ,rcomhl]=Fresnell(f ,psi) ; 

EPSREL=8;SIGMA=.04; 

[rcomv2 , rcomh2] =Fresnell (f ,psi) ; 

EPSREL=80;SIGMA=.004; 

[rcomv3 , rcomh3] =Fresnell (f , psi ) ; 

EPSREL=2 ; SIGMA= . 000 ; 

[rcomv4 , rcomh4] =Fresnell (f , psi ) ; 

subp lot ( 2 , 2 , 1 ) , plot ( ( 180/pi ) *psi , abs (r  comv) ,  ’  k  ’ ) 
hold 

axis( [0,90,0,1] ) 

xlabel( *  Incidence  angle  (degrees) O 
ylabel( ^Reflection  Coefficient’) 

'/,text(5, .8, [’SIGMA  =  ’ )num2str (SIGMA), ’  S/m  EPSREL=  ’,num2str(EPSREL)]) 

subplot ( 2 , 2 , 1 ), plot (( 180/pi ) *psi , abs (rcomvl ),’ r  ’ ) 

subplot (2,2,1), plot ( ( 180/pi ) *psi , abs (r comv2) , 1 b J ) 

subplot(2,2,l),plot((180/pi)^psi,abs(rcomv3),’g’) 

subplot (2,2,1), plot (( 180/pi ) *psi , abs (r comv4) , J  c J ) 

hold 

subplot(2,2,2),plot((180/pi)*psi,(180/pi)*angle(rcomv), ’k’) 
hold 

axis(  [0,90,-190,5] ) 

xlabel( ’Incidence  angle  (degrees)’) 

ylabel(’Phase’) 

subplot(2,2,2) ,plot( (180/pi)*psi , (180/pi)*angle(rcomvl) , ’r’) 
subplot (2, 2 ,2), plot((180/pi)*psi, (180/pi) *angle(rcomv2),,b’) 
subplot(2,2,2) ,plot( (180/pi)*psi , (180/pi)*angle(rcomv3) , ’g’) 
subplot(2,2,2) ,plot( (180/pi)*psi ,-(180/pi)*angle(rcomv4) , *0') 
hold 


subplot(2,2,3) ,plot((180/pi)*psi ,abs(rcomh) , ’k’) 


hold 

axis( [0,90, .5,1] ) 

xlabeK1  Incidence  angle  (degrees)1) 
ylabel( 'Reflection  Coefficient’) 
subplot(2,2,3)  ,plot((180/pi)+psi,abs(rcomhl) ,  "rO 
subplot(2,2>3) ,plot((180/pi)*psi,abs(rcorah2) , 'b1) 
subplot (2 , 2,3) , plot ( ( 180/pi ) *psi , abs (rcomh3) ,  ’  g  ’ ) 
subplot(2,2,3) ,plot((180/pi)*psi ,abs(rcorah4) , ’c') 
hold 

subplot(2}2,4)  ,plot((180/pi)*psi ,  (180/pi)+angle(rcorah) ,  T ) 
hold 

axis(  [0,90,170,200] ) 

xlabeK’ Incidence  angle  (degrees)’) 

ylabel( ' Phase’) 

subplot(2 ,2,4) ,plot((180/pi)*psi , (180/pi)*angle(rcomhl) , ^r’) 
subplot(2,2,4) ,plot( (180/pi)*psi , (180/pi)*angle(rcomh2) , ’bO 
subplot(2,2,4),plot((180/pi)*psi, (180/pi)*angle(rcomh3) , Jg’) 
subplot(2,2,4)  ,plot(  (180./pi)*psi ,  (180/pi)*angle(rcomh4) ,  ’c J) 


hold 
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